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Abstract. Flat directions in generic supersymmetric theories can change the thermal 
history of the Universe. A novel scenario was proposed earlier where the vacuum 
expectation value of the flat directions induces large masses for all the gauge bosons 
and gauginos. This delays the thermalization of the Universe after inflation and solves 
the gravitino problem. In this article we perform a detailed calculation of the above 
scenario including certain aspects not considered earlier and show that the gravitino 
problem is actually not solved for certain initial values of the flat direction vacuum 
expectation value. In fact the large gaugino mass can lead to enhanced gravitino 
production. 
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1 Introduction 

Large vacuum expectation values (vevs) for fiat directions in supersymmetric (SUSY) 
theories, i.e., directions in the field space of scalars that have a flat potential, can affect 
the thermal history of the Universe in two ways as discussed in the literature. Firstly, 
they may delay the thermalization of the inflaton decay products [1-3]. Secondly, 
the flat directions could potentially dominate over the inflaton decay products [4, 5], 
such that reheating is associated with the decay of flat directions as opposed to the 
decay of the inflaton. Below we shall consider gravitino production in the context of 
flat directions delaying thermalization of the decay products of the inflaton which are 
presumed to dominate the Universe. 

Gauge symmetries are broken along flat directions and for certain flat directions 
all standard model (SM) gauge symmetries are broken in the early Universe. (Low 
energy phenomenology is not affected as the vevs of these flat directions ultimately 
disappear.) In refs. [1, 2] the authors have presented an interesting scenario in which 
flat directions give a SUSY conserving large mass ~ yfaip to all gauge bosons and 
gauginos where tp is the vev of the field <j) that parametrizes the flat direction. Since 
processes involving gauge bosons are important for thermalization, this results in a 
slow rate of thermalization after the end of inflation. Thus the Universe enters a 
period of pre-thermalization after inflation. Throughout this article we use the term 
thermalization to refer to chemical equilibration, and thermal equilibrium to refer to 
chemical equilibrium. 

In the scenario of refs. [1, 2] after inflaton decay the Universe reaches a state of 
kinetic equilibrium via 2—7-2 interactions after a certain time interval has elapsed. 
Further thermal equilibrium is reached only when number violating interactions, i.e. 
2—7-3 interactions, are effective. Thus the rate of inflaton decay, of kinetic equilibration 
and of thermal equilibration follow the relation 



T d 3> Tk in > ^thr i 




where 



T d = Inflaton decay rate 
T^in — Rate of kinetic equilibration 
T t hr = Rate of thermal equilibration, 
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and we define tj, tk in and t t hr as the times of inflaton decay, kinetic equilibration and 
chemical equilibration respectively. 

In refs. [1, 2] the authors argue that due to the delay in thermalization the generic 
problem of gravitino overproduction in SUSY models is also avoided. The lack of 
thermalization leads to a dilute plasma of high energy particles, and the low number 
density of particles effectively leads to low gravitino production. Moreover due to a low 
final reheat temperature after thermalization gravitino production after thermalization 
is also suppressed. 

To find the number density of gravitinos ng, refs. [1, 2] use the integrated Boltz- 
mann equation of the form 



from inflaton decay till thermal equilibration for the process / + J — > G + K, where / 
and J are quarks or squarks and K is a gauge boson or a gaugino. We too obtain and 
use this equation for t kin < t < t thr , where one can describe the particle distribution 
function by a thermal distribution with an effective chemical potential. For t < tj, in 
we explicitly evaluate the collision integral in the integrated Boltzman equation using 
a delta function distribution function for the incoming particles. 

There are two important effects associated with large vevs of the flat directions 
as relevant for gravitino production that have not been considered earlier. 

1. Gauge bosons and gauginos are produced along with gravitinos in the relevant 
processes for gravitino production. As mentioned above these get large masses 
due to the large vevs of the flat directions. Consequently gravitino production 
ceases in cases where the energy of quarks and squarks in the Universe falls below 
the mass of gauge bosons and gauginos due to phase space suppression. 

2. The large m~ g increases \M\ 2 in the cross-section for gravitino production. This 
can greatly enhance the gravitino abundance. 

As we show below, both these effects can have important consequences for the total 
gravitino abundance. (In ref. [6] we have included these effects in a study of gravtino 
production in a thermal Universe.) As in ref. [2] we consider two cases, mo < and 
tuq > Td, where mo is the condensate mass. We find that for mo < phase space 
suppression shuts off gravitino production both before thermalization and even after 
thermalization till the condensate decays. For m > the gravitino abundance is 
greatly enhanced due to the the large gluino mass and is several orders of magnitude 
larger than the cosmological bounds on the gravitino abundance. 

It should be noted here that the longevity of these flat directions has been in 
dispute [3-5, 7-13]. One of the key issues is whether or not non-perturbative effects lead 
to a fast decay of the condensate. In ref. [7, 8] it was pointed out that non-perturbative 
decay due to parametric resonance is suppressed for a complex scalar condensate if the 
real and imaginary parts of the scalar field oscillate out of phase. In that case the 
field does not pass through the minimum of the potential at ip = and so the required 




(1.3) 



i<j 
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non-adiabatic condition \ujk\ ~ w\ is not satisfied, where ojk is the the frequency of 
the /c-th mode of the x field produced by the decay of the condensate. (Backreaction 
effects can also be important in suppressing non-perturbative decay.) The condensate 
then decays only perturbatively and hence lasts for a long time. Ref. [4] however 
argued that in the case of multiple flat directions there is a mixing of excitations 
along different directions which makes non-perturbative effects important. Further 
investigations in ref. [9] implied that even with multiple flat directions resonant decay 
is most likely to be suppressed (in cases of physical interest) because the multiple flat 
directions effectively reduce to one flat direction or a collection of independent single 
flat directions. Raising certain issues related to gauge fixing, refs. [10, 11] worked 
in the unitary gauge to eliminate the Nambu-Goldstone boson that appears in the 
context of SUSY flat directions charged under the gauge group of the MSSM and 
studied the circumstances under which the presence of multiple flat directions can lead 
to non-perturbative decay. Ref. [12] invoked charge conservation to argue that non- 
perturbative decay is suppressed. At the same time ref. [5] analysed the issue with 
proper gauge fixing and concluded that with one flat direction there is no resonant 
particle production, in agreement with ref. [12], while disagreeing with the conclusion 
of ref. [12] in the context of multiple flat directions. However refs. [3, 12] argue that even 
if non-perturbative decay occurs for multiple flat directions it leads to redistribution 
of energy of the condensate amongst the fields in the D flat superspace and hence 
to practically the same cosmological consequences (including large gauge boson and 
gaugino masses and delayed thermalization) as in the scenario with only perturbative 
decay. 

Our analysis below presumes that perturbative decay is the relevant decay mech- 
anism of the flat direction condensate and that the impact of the condensate lasts long 
enough to have cosmological consequences. This may be because either there is only a 
single flat direction that is excited, or because the multiple flat directions are excited 
such that parametric resonance is still suppressed, or because of the arguments of refs. 
[3, 12], or because of backreaction effects. 

There are other effects that affect the evolution and decay of the flat direction 
condensate. Thermal contributions to the condensate potential can change the mass 
of the condensate leading to an earlier oscillation time for the condensate [14-16]. 
Below we take the condensate mass mo = 100 GeV . We present our results in terms of 
mo which can be replaced by larger values to include these effects. (Ref. [16] obtains 
masses of ~ 10 10 GeV , 10 6 GeV and 10 5 GeV for n = 1, 2, 3 where terms that lift the 
flat directions have the form 2n+4 /M^ n where M* = 10 18 GeV .) Another effect is that 
the condensate can decay due to scattering by inflaton decay products [14, 15, 17]. 
Both these effects have been discussed in appendix A. 4 of ref. [2] where it has been 
argued that these effects will not be relevant during the non-thermal phase. 

The flat direction condensate can also decay to Q-balls due to inhomogeneities in 
the condensate field. This has been studied in the context of both gauge mediated and 
gravity mediated SUSY breaking [18-26]. In addition it has been seen that condensate 
lumps with different properties than Q-balls called Q-axitons may also be formed [27]. 
Note however that the time scale for the formation of Q-balls and Q-axitons can be 
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large. For gauge mediated SUSY breaking scenarios the time scale for Q-ball formation 
t ~ 5 x 10 5 m _1 [20], where m is a characteristic mass scale for the flat direction 
condensate, and for gravity mediated SUSY breaking scenarios r ~ 100 — 1000 m~ l 
[19] or 5.5 x 10 3 m _1 [21]. Ref. [18] also obtains a Q-ball formation time which is larger 
than the flat direction oscillation time by a factor of 1.5 x 10 4 . The time scale for 
Q-axiton formation is 1600 m -1 [27]. Q-ball formation may be avoided for certain flat 
directions, including those with a large stop admixture, which allow for K > 0, where 
K is a parameter in the potential in gauge mediated SUSY breaking scenarios [22]. 
Our analysis does not include Q-ball or Q-axiton formation. 

We now present a summary of this article. In section 2 we provide the integrated 
Boltzmann equation and phase space distribution functions valid for a non-thermal 
Universe. Section 3 provides the scattering processes relevant for gravitino produc- 
tion in the scenario with heavy gauge bosons and gauginos. In section 4 we obtain 
the relevant time scales, namely, tki n , t t hr, tf and to where tf is the time of decay of 
the flat direction condensate and to is the time in the radiation dominated era after 
thermalization when gravitino production could commence. In this section we also 
discuss the feasibility of gravitino production in the context of phase space suppression 
in different epochs. The cross-section for gravitino production in the center-of-mass 
frame has been calculated in section 5. Then in section 6 we obtain the gravitino 
abundance generated in different time intervals in the Universe. The parameters re- 
quired to obtain the gravitino abundance are given in section 7. Our results are given 
in section 8 and our conclusions in section 9. The Appendices contain derivations of 
certain expressions used in the text. 

2 The Boltzmann equation 

To obtain the number density of a species X 3 participating in reactions X\X 2 ^ X 3 X 4 
the integrated Boltzmann equation is 



gi is the number of internal degrees of freedom of species i. The |.M| 2 includes an 
average over initial and final internal degrees of freedom [28]. 

Typically X\ and X 2 are assumed to be in full thermal equilibrium with zero 
chemical potential and thus the phase space distribution function has the form 





where fa are phase space distribution functions and 



Hi d 3 pi 



(2.2) 





i = 



1,2. 



(2.3) 
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This is because X\ and X 2 usually have other stronger interactions than their interac- 
tion with Xs, X4. This assumption is not valid in the scenario considered in refs. [1, 2] 
as the Universe is not in thermal equilibrium before tthr- bn fact /i,2 have the following 
form 

f i = C i 6(E i -'^-^) for t d <t<t km (2.4a) 



2 a 
E - & 

exp ( — - ) for t kin <t < tthr (2.4b) 



where, 



Ci =e*&(^)f®-, i = l,2 (2.5) 



ml \a / 7rp. 



as derived in appendix A. In the above expressions is the mass of the inflaton, 
pcftitd) is the energy density of inflaton at time t d , a is the scale factor and is 
the effective chemical potential. Eq. (2.4a) reflects the scenario that before tkin the 
relativistic energy of fermions and sfermions (whose interactions produce gravitinos) is 
their energy at the time of their production from (2 body) inflaton decay at td scaled 
by 1/a. Between t^m and t t h r one does not have full thermalization as represented by 
& in eq. (2.4b). (Also see footnote 24 of ref. [2].) 

When the number density of X 3 is small, as in our case where X 3 represents the 
gravitino, we can ignore the product f 3 f± in eq. (2.1). Thus, eq. (2.1) can be written 
as 

h 3 + 3Hn 3 = J tflli dU 2 h h = A, (2.6) 

where we have defined a variable A to denote the collision integral on the right-hand- 
side of the above equation. Here, W\2{s) is dimensionless and Lorentz invariant, and 
dependent only on the Mandelstam variable s. W12 is given by 

W 12 (s) = j dTl 3 dTl 4 (2ir) 4 5 4 (p 1 +p 2 -p 3 -p 4 )\M\ 2 } (2.7) 

which can also be written as AE\E 2 ov re \ where v re i is [(pi-Pi) 2 — mlm 2 ,] 1 ^ 2 /{E\E%). 
2.1 The collision integral 

For the completeness of the paper we provide below some details on how to obtain the 
collision integral. It is easier to calculate W\ 2 of eq. (2.6) in the centre-of-mass (CM) 
frame and it can then be expressed as [29] 

W l2 {s) = Ap 12 a CM {s) , (2.8) 

where 

[s - (mi + m 2 ) 2 ] 1/2 [s - (mi - m 2 ) 2 ] 1/2 

P12 = ( 2 - 9 ) 
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is the momentum of particle X\ (or X 2 ) in the center-of-mass frame of the particle 
pair (Xi,X 2 ). The cross-section in the CM frame oqm is 



<?cm{s) 



dQ 



\dn) 



(2.10) 



CM 



where 



da 

dn 



CM 



9394 
6AtT 2 S 



[s - (m 3 + m 4 ) 2 )(s - (m 3 - m 4 ) 2 
_(s - (mi + m 2 ) 2 )(s - (mi - m 2 ) 2 )_ 



1/2 



I-MP 



(2.11) 



is the differential scattering cross-section in the CM frame [30]. Note that the averaging 
over the final state degrees of freedom in \Ai\ 2 cancels the prefactor g 3 g 4 in the above 
equation. \A4\ 2 (s,t,u) in eq. (2.11) can expressed in terms of s and 6 13 , the angle 
between X\ and X 3 in the CM frame, using u — —s — t + £\ m 2 and t - 
2E X E 3 + 2a/-E 2 - m\ y 1 E\ - m\ cos9 13 , with E x = \Jp 2 12 + m\ and E 3 
ml)/(2 v rs) [30]. 

To perform the integration over dU\ dn.2 in (2.6) we note that 



m\ + m 2 — 
= (s + m 2 . — 



dlli dn.2 = 



91 



d 3 pi 



(2?r) 3 2^1 

and that the volume element d 3 pi d 3 p 2 is given by 



92 



d 3 p$ 
(2vr) 3 2^ 



(2.12) 



d 3 p 1 d 3 p 2 = 4:7i\p 1 \E 1 dE 1 4:TT\p 2 \E 2 dE 2 - d(cos 9 21 ] 



(2.13) 



where #21 is the angle between X\ and X 2 . (Note that we now work in the rest frame 
of the Universe.) In order to simplify the integration we replace the variable cos 6 2 i 
by s via the relation s = m\ + m\ + 2EiE 2 — 2\pi\\p 2 \ cos #21- The lower limit of the 
integral over s should not be less than either (mi 
the phase space factor in eq. (2.11). Therefore 



m.2) 2 or (m 3 + m 4 ) 2 as dictated by 



dlli dU 2 



9i 92 
4(2tiY 



dEx 



nil 



dE 



111-2 



ds, 



(2.14) 



where s m i n = max [s±, (mi + m 2 ) 2 , (m 3 -|-m 4 ) 2 ] and si^ 



m; 



■m 2 2 + 2E 1 E 2 T2\p l \\p 2 \. 



For mi,m 2 ,m 3 « 0, s min m m 2 . Using eq. (2.14) in eq. (2.6) we get 

9i 92 



h 3 + 3Hn 3 



4(2tt) z 



dE 1 



nil 



dE, 



in 2 



S2 



dshf 2 W l2 . 



(2.15) 



Thus, after calculating \AA\ 2 for each process mentioned in the next section, we 
follow the procedure given above to find the parameter W\ 2 (s) and then with the fi 
as given in eq. (2.4) we use the above equation to calculate the number density of 
gravitinos produced between t d and tki n , and t kin and t t h r - 
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3 Relevant processes for gravitino production 



The production of gravitinos after inflation has been discussed in [31-50] The processes 
involving strong interactions denoted by A to J are [35] 

• A: g A + g B ->■ g c + G 

• B: g A + g B -> g c + G 

• C: ft + g A ->■ ft +G,^ + g A ->■ ft + G 

• D : q% J rg A - J f ft + G, ft ± g A ->■ <| + G 

• E: ft + |j+ -> £ A + G, ft + ft ->■ # A + G 

• F: ~g A + # B -> ^ c + G 

• G: ft + g A -> + G, ft + ^ A — )■ ft + G 

• H: ft + <? A — > ft + G, ft + g A -> ft + G 

• I: ft + ft — )■ ^ A + G 

• J: ft + ft ->• ^ A + G 

The above can be generalized for all gauge groups. We shall consider SU(3) C only. In 
the non-thermal scenario we only have to consider processes involving the scattering of 
fermions and/or sfermions [2]. The reason is as follows: the gauge bosons and gauginos 
have large masses induced by the flat direction vev. Due to their large masses they 
decay quickly to lighter fermions and sfermions and therefore do not participate in the 
production of gravitinos. The processes not involving gauge bosons and gauginos as 
incoming particles are the E, I and J processes. These processes are shown in figures 
1, 2 and 3. (The second E process is not shown and its contribution will be the same 
as that of the first process.) 

In ref. [42] it has been argued that the gravitino production cross-section should 
be proportional to 

(31) 

where M = Mpi/\/8tt = 2.4 x 10 18 GeV is the reduced Planck mass. The factor 1 is 
related to the contribution of the ±3/2 helicity states of the gravitino while the term 
proportional to m|/(m~) represents the contribution of the helicity ±1/2 states. To 
obtain the contribution of the ±1/2 helicity states, for s/s ^> mg we may express the 
gravitino field as 

^Vfi 8 "* (3 - 2) 

where ip represents a spin 1/2 fermionic field. Cross-sections for spin 1/2 gravitino 
production can then be obtained from the corresponding effective lagrangian for ijj 
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Figure 1. Tree level diagram for E: qi + qj — >■ g A + G scattering. The figures for + q~j — > 
g + G are similar. 




(a) (b) (c) 

Figure 2. Tree level diagram for I: qi + qj g A + G scattering. 

[51-53]. The total cross-section is then obtained in ref. [42] using eq. (3.1). In the 
scenario under consideration, m g ^> mg and the contribution of the ±1/2 helicity 
states dominates. 
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Figure 3. Tree level diagram for J: + qj — > g A + G scattering. 

4 Obtaining tunt Uhn an d tj and checking the feasibility of 
gravitino production 

In this section we obtain tu n an d tthr for the cases m < Td and mo > Td considered in 
ref. [2]. Since gravitino production is suppressed by lack of thermalization we impose 
the constraint that the Tkin,thr < H at to study this effect. We also impose the 
constraint that the energy of the incoming quarks and squarks should be greater than 
3m Si j to avoid phase space suppression of the gravitino production cross-section by 
the s — rn 2 g g factor in eq. (2.11). These constraints give lower and upper bounds 
on ip respectively. In the case of m < Td, these bounds are in fact incompatible 
indicating that lack of thermalization effectively shuts off gravitino production through 
phase space suppression. Surprisingly, we find that even after thermalization, for both 
cases mo ^ Td, gravitino production is phase space suppressed till the flat direction 
condensate decays because of the large gluon and gluino mass. 

We take the inflaton mass m^ to be 10 13 GeV and a = gf/(47r) to be 5 x 10~ 2 
as applicable for the relevant energy scales for t < t t h r and a = 0.1 for t 3> tthr as 
relevant below. We presume the inflaton decays via (j) — >■ X\ + X 2 so that the number 
density of the X\ or X 2 particles equals that of the inflaton. The time when the flat 
direction starts oscillating is to = m o \ an d the inflaton decays at td = T^ 1 . The gluon 
and gluino mass is given by 

2 J a( P0 fOT * < *0 /a -j \ 



a(fiQ (do /a) for t > t 



o 

4.1 m < T d 

For mo < Td, the flat direction condensate starts oscillating after the inflaton decays. 
For this case we take m = 100 GeV and Td = 10 4 GeV as in section 3.1 of ref. [2]. 
One must separately consider the time intervals td — t , t —t kin , tk in — t thT and t > t thr . 

The gluon and gluino masses are given by \fatp. The energy of the incoming 
particles is m^dd/a) for td < t < tkin- Since kinetic equilibration does not change 
the number of particles while redistributing their energies, the average energy of the 
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incoming particles for t k i n < t < tthr is also given by m^dd/a). Gravitino production 
for t > tthr is considered in section 6.1. 

4.1.1 Non-equilibrium condition 

At t = td we impose 

H{t d ) = T d > r kin (td) = n(av) kin = ^Ml!^! (4.2) 

where we have used p^itd) = (3/(87r))Mp^ Y 2 d and {av) k i n ~ 4:iia 2 /m 2 [1, 2]. Then we 
get 

3M|T> 

J- d > 



2 ifl 
3fTd^ l/2 



or, ^>y-^_J 
or, <^ >10 14 GeV, (4.3) 

where we have used m 2 = onp\. H(td) > F thr (t d ) gives ip > 2 x 10 13 GeV. 
Note that for t d < t < t , T kin is given by 

Tfcm = n(av) kin 

pSd) fa d \ 3 4:7ia 2 ^ 

(4.4b) 



1 1 

a 3 t 3 ' 2 ' 



where we have used n = n(td){a,d/a) 3 ■ Since H ~ we see that is falling faster 
than H. So, if H > T kin holds at td, then it also holds for td < t < t . 

4.1.2 Equilibration times 

For t > t , T kin is given by eq. (4.4a) with m 2 = a(pl(ao/a) 3 . Therefore 



3 M 2 Pl T 2 /aA 3 a_ 
2 m \a Q ) <Po 



Tu^V-^r-) ~ 2 (4.5a) 



= ^77 rV 2 m^ 2 a = constant. (4.5b) 

Since if is decreasing with time, after some time t k i n > to, H becomes less than T k i n , 
and later it becomes less than T t hr- Setting t kin = T^ n we get 

or, t kin = 9 x 1(T 30 GeV - 1 • (4.6b) 
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rW(t) = n(av} thr = aT kin as (av) thr ~ 47ra 3 /m 2 [1, 2]. For t > £ , F thr (t) ~ 
(1/a 3 ) (a 3 /m 2 ) ~ constant, as for T kin . Then t t hr = a~ l t kin . 

As we shall see below, there is also an upper bound on tpo which is highly restric- 
tive. Taking (p = 10 14 GeV , we get, 

t kin = 9x lO^GeV- 1 . (4.7) 

Note that t d = lO^GeV" 1 and t = 10 -3 GeV -1 . t thr = 2 GeV" 1 . For larger <p , 
both t kin and tthr will be larger. 

4.1.3 Avoiding phase space suppression 

Gravitinos are produced along with gluons or gluinos in the processes discussed in 
section 3. To avoid phase space suppression by the s — m 2 ^ factor in the differential 
cross-section we impose 

3m sj < E x + E 2 . (4.8) 

We apply this condition for gravitino production before and after thermalization. 

Before thermalization: We first consider gravitino production in the time interval 
t d <t <t . Eq. (4.8) implies 

< (—) . (4.9) 



a 

where E\ and E 2 are the incoming energies in the CM frame. We see that the RHS 
falls like 1/a, whereas the LHS is constant. We wish to analyse whether delayed ther- 
malization suppresses gravitino production. Therefore to consider maximum gravitino 
production, we impose the above inequality at to- If the above inequality holds at to, 
then it will hold at earlier times till t d too. Using the above inequality we find an 
upper bound on <p . 

1/2 



or, 



rn^ ( t d 
3y/a \t 

^° \T d J 



1/2 



or, (f < 1.5 x 10" GeV. (4.10) 

Thus we need cpo < 1.5 x 10 12 GeV to get gravitino production from t d till time to- For 
the time interval t < t < t t h r , m g ^ ~ ~ 1/a 3 / 2 and the condition in eq. (4.8) gives 

3^ (^) 3/2 <^(^). (4.11) 

We see that the RHS falls like 1/a and the LHS falls like 1/a 3 / 2 . So if the above 
inequality holds at £ 0; which implies eq. (4.10) then it will hold for all time t till tthr- 
We see that for the time interval t d < t < t thr the upper bound on (p (eq. (4.10)) 
is smaller than the lower bound on <p (eq. (4.3)). Thus, gravitino production for the 
entire time interval t d < t < tthr is n ot possible for the mo < T d case. 
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After thermalization: For t > tthr-, i-e. after full equilibration, the average energy 
of particles in the Universe falls. Gravitino production starts at to when 3m g ^ becomes 
less than the incoming particle energy (E\ + E 2 ), i.e., 



3m ffi j < E 1 + E 2 



or, 



(4.12) 



where the subscript R denotes reheating or full equilibration at t t h r - This implies 



(—) ( 



3/2 fa thr \ 3/2 



or, 



Hhr 



or, 



JH>0 Uhr) 



-3/2 



Then, 



3/4 



a thr 

a 
a 

O'thr 



< 



ttthr 



1/2 



(4.13) 



n 



j_\ 3/2 

2,„4 r2 



81 (it) 2^21 



3 ' 
'"O 



(4.14) 



where we have used section 4.1.2 for tthr and the expression for Tr is derived later in 
section 7.5. A is obtained to be 4 x 10~ 5 in section 7.3. 

The perturbative flat direction decay rate is ~ m\j^) 2 [4, 54] and the flat direction 
decays when the decay rate is of order H at a time tf. Then 



or, 



t _ 
J / ~~ — 3 
rriQ 

= ^0 
mo 

^0 

4/5 
9/5 



(4.15) 



a f 



3/2 



(4.16) 



rn 



We see that to < tf only for <£>o < 3 x 10 12 GeV, which conflicts with the lower 
bound on tp in eq. (4.3). Thus the non-equilibrium condition implies that the conden- 
sate decays before gravitino production can commence in the radiation dominated era 
after thermalization. Subsequent to the decay of the condensate there will be gravtino 
production in the thermal Universe as discussed later in section 6.1. 
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4.1.4 Relaxing the bound from phase space suppression 

If we relax the bound in eq. (4.10), we may get gravitino production for some part of 
the time interval t d < t < t thr . We first consider < t < t . Using a ~ t 1//2 in eq. 
(4.9) and td = T^ 1 we find gravitino production is not phase space suppresed for 

t m l 1 
T d 9a<f 2 

or, t < 2 x 10- 6 GeV _1 , (4.17) 

for ipo = 10 14 GeV. Since this corresponds to times less than td = 10~ 4 GeV _1 , it 
implies that one cannot get gravitino production anywhere in this time interval. 

We now consider t < t < tthr- Gravitino production can start at some time after 
to, when the inequality in eq. (4.11) is satisfied. This is equivalent to 

or, Sla'^f^-X <t. 
m m^ \m J 

For ip = 10 14 GeV this implies 

t > 2 x 10 5 GeV -1 . (4.19) 

Since this is larger than t t h r = 2 GeV -1 , once again there is no gravitino production 
in this time interval. 

The above discussion shows that the conditions required for lack of kinetic and 
chemical equilibration leads to a large outgoing gluon and gluino mass and thus to phase 
space suppression that shuts off gravitino production for the entire time interval from 
td to t t h r in the m < Td case. Increasing the value of (fo will not affect this conclusion 
as this will only lower or raise the limits in eqs. (4.17) and (4.19) respectively. Similarly, 
gravitino production is phase space suppressed in the radiation dominated era before 
the condensate decays at t/ and increasing the value of <po will not alter this conclusion. 
(Decreasing the value of y?o to avoid phase space suppression will lead to the standard 
scenario of reheating and gravitino production in a thermal Universe.) 

Thus, for the parameters chosen above for the scenario under consideration, there 
is no gravitino production for the case m < Td before the condensate decays. 

4.2 m > T d 

We now consider the case m > Td, where the flat direction condensate starts oscillating 
before the inflaton decays. m ~ 100 GeV and we take T d = 10 GeV as in section 3.1 
of [2]. 



(4.18) 



4.2.1 Non-equilibrium condition and relevant time scales 

The absence of kinetic equilibration of the inflaton decay products requires Tk in < 
H. Tkin is given by eq. (4.5a). We require Tkin to be much less than H at t d to 
obtain a significant phase of non-equilibrium since H decreases with time while Tkin 
remains constant. Then applying the non-equilibrium condition at td, and setting 
= vl^l/al) and H{t d ) = T d , we get 

3 ( ^ ^ 

or, tp > 3 x 10 13 GeV , (4.20) 

where we have used a ~ t 2 ^ 3 between to and td as relevant for the inflaton oscillating 
in a quadratic potential before decay. Chemical non-equilibration at td requires </?o > 
3 x 10 13 GeV ja~7x 10 12 GeV . 

Thus, for (po < 7 x 10 12 GeV the Universe is in kinetic and chemical equilibrium 
at t d . For 7 x 10 12 GeV < y? < 3 x 10 13 GeV the Universe is in kinetic but not in 
chemical equilibrium at td. Finally, for <£>o > 3 x 10 13 GeV the Universe is neither in 
kinetic or chemical equilibrium at td- 

To obtain t k i n we use eq. (4.5a) and a ~ t 2 / 3 for t < t < t d . Then 

! 3 aM 2 ,r 2 (t d 

l kin — L km — 9 2 1 + 

I m^ipQ \t 
3aMj, t T 2 d fm 



m 4> <Po V r 



d 



3 a M 2 Pl m 2 



2 m^ipl 

As before, r tfer = aT k i n and i thr = T^ r = a~H kin . 



(4.21) 



4.2.2 Avoiding phase space suppression 

For gravitino production not to be phase space suppressed by the s — m 2 g factor for 
td < t < t t hr we again impose the condition of eq. (4.8) and get eq. (4.11). To consider 
maximum gravitino production we assume that this condition is valid from td- Using 
eq. (4.11) at td gives 

3Vatp d < m 4 , (4.22) 
where tp d = <po{a Q /a d Y /2 = tp (t /t d ) = Lp (T d /m ). Therefore we get 

<p < 1.5 x 10 14 GeV. (4.23) 
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Since the LHS of eq. (4.11) falls faster than the RHS, this upper bound on ifo also 
ensures the viability of gravitino production at times after td. 

We shall take <po = 10 14 GeV hereafter so that gravitino production commences 
at td and there is neither kinetic nor chemical equilibration at t d - For this value of ipo, 
t kin = 0.9 GeV - 1 and t thr = 20 GeV - 1 . 

For gravitino production after tt hr , eq. (4.12) implies 



ad J 



3/2 



Od_ 

Clthr 



3/2 



0>thr 



3/2 3 

<2x-T R 



or, 



V®<Po 



or, 



<*<Po 



m 



td_\ 

Uhr J 



3/4 



a 



0>thr 



0>thr 

a 

1/2 



3/2 



r<i tthr J 



< 



dthr 



(4.24) 



Then 



tt (s: 



,1/2' 



rrir 



1 

tthr 



2,^4 F3 
<2 



(4.25) 



where we have used the expression for tthr from section 4.2.1 and for Tr as in eq. (7.8) 
or (7.10). A = 4 x 10" 11 as obtained in section 7.3. 

The time when the condensate decays, tf is given by eq. (4.15). Then 



or, 



*/=4i 

mf. 



t 



f 



^0 

m„ 



a d 

to 
t d 

m 

4/5^1/5 



3/2 



3/2 



(4.26) 



For <p = 10 14 GeV, t G is 5 x 10 11 GeV" 1 while t f is 3 x H^GeV -1 . Thus, there 
is no gravitino production in the radiation dominated era after thermalization before 
the flat direction condensate decays. 



4.2.3 Relaxing the bounds 

The lower and upper bounds on ip in eqs. (4.20) and (4.23) are restrictive. For later 
use we relax the upper bound from phase space suppresion. If <po > 1.5 X 10 14 GeV it 
implies that gravitino production commences at a time t p > td. To get an expression 
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for t p we start with eq. (4.11) and solve for a p , the scale factor at which gravitino 
production starts. We thus get 




(4.27) 



where we have used the fact that for t < t < td the scale factor a ~ t 2 / 3 and for 
td < t < t p , a ~ t 1 / 2 . One would then get an extra factor of 1/(1^ t p ) 7 / 2 along with td 
replaced by t p in the expressions for gravitino number density generated before kinetic 
equilibration in eqs. (6.5a), (6.8a) and (6.10a). 

We remark in passing that if 7x 10 12 GeV < ipo < 3x 10 13 GeV and one has kinetic 
equilibration of the inflaton decay products at td without chemical equilibration, then 
one would obtain gravitino production using only section 6.2.2, with T kin set equal to 



5 Calculations of a cm 

Having discussed the feasibility of gravitino production in a non-thermal Universe in 
the context of phase space availability we now calculate the cross-section of gravitino 
production. The expressions for \Ai\ 2 in this section include an average over spins and 
colors of the initial and final state particles but are for a specific flavor for quarks and 
squarks. 



5.1 Differential scattering cross-sections 

The differential scattering cross-sections in the CM frame is as given in eq. (2.11). We 
obtain \M\ 2 with help of MATHEMATICA package FeynCalc. We assume that the 
mass of the gravitino is small compared to the mass of the inflaton and approximate 
it to be zero in the numerator of a cm- Also the left-right squark mixing angle 9 is set 
to zero. The \M\ 2 for helicity ±1/2 contribution for the E, I and J processes are given 
by 



\Mi 



V n 2 \T A \ 2 

/L^A,i,j vs I x 3% I 

3#i 92 9s 94 M 2 s 4 m| 



2s(m 2 + s) 2 m 6 g + 2(m 2 g + s) (t 



m 



4 3s 2 m 2 



+ s (-At ml + 2s(3s - 2t) m\ + 4s 3 m~ + sH) 
+ 2s 2 m 2 s (tm~ + s(t - s) m~ + s 2 t) 



(5.1a) 



\Mi 



j 1 1 



3#i 92 93 giM 2 s 3 m| 



(( m 



s) (m| - 2(s + t) m 2 + s 2 + 2t 2 + 2st) 



+ 2s m 2 g (ml - 2(s + t) m| + s 2 + 2t 2 + 2st) 
+ s 2 m? ( m| - (s + 2t) m\ + 2i(s + t)) 



m. 



(5.1b) 



J* I 



3^i 52 53 g4M 2 s 3 m| 



+ 2s m~ m 2 + s 2 m~ 



( m? - s) s 3 + 2t ( m 2 - s - t) ( ( m\ - s) 



m. 



(5.1c) 



where T A is the generator of the SU(N) group in the fundamental representation, and 
i, j represent the color of the incoming particles. If we consider only terms quadratic 
in m g and rrig then we get 



\Mi 



\Mj\ 2 



l-M 



9 2 s E^m||^| 2 i-t)(m 



2m?) 



M 2 3g 1 g 2 g 3 g 4 m 2 s 
M 2 ?>g x g 2 g 3 g±m 2 ~ 



It) (s + t) 



g 2 E^jg^T 2(s 2 + 2ts + 2t< 
M 2 3 g x g 2 g 3 #4 m% s 



(5.2a) 
(5.2b) 
(5.2c) 



We see that [.A/f^jl 2 agree with the helicity ±1/2 contribution of Table 1 of ref. [42], 
while |A^_b| 2 has a contribution from the massive gluon. 1 

We shall work with the full amplitude squared eq. (5.1) and only choose the 
dominant terms in the limit m^dd/a) ^> mg, m g ~ g when we consider the A parameter. 
We shall take m 2 = m~ because the mass due to the SUSY flat direction is SUSY 
conserving [2]. 

The differential scattering cross-sections for the three processes are 



2s ( m~ + sf m 6 g + 2 ( m~ + s) (t m\ - 3s 2 m~ - s 2 {s + t)) m\ 

+ s (-At m~ + 2s(3s - 2t) m| + 4s 3 m\ + s 3 t) m 2 + 2s 2 m\ (t m\ + s(t - s)m 2 g + s 2 t) , 

(5.3a) 

1 |7\4/| 2 above differs from ref. [42] by a factor of 2 because we sum over both chiralities of the 
quarks. The expressions in Table 1 of ref. [42] do not include an average over incoming and final 
states and hence we have a sum over A, i,j and an additional factor of g± gi g% <?4 in the denominators 
of \M e ,i,j\ 2 . 
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daj v,,,,//;('S' *) \t- 



A\2 



48 g 2 M 2 vr 2 s 5 



HI} 



s) (m 



2(s + t) m 2 s + s 2 + 2t 2 + 2st) 



m„ 



+ 2s m\ ( mj - 2(s + t) mf + s 2 + 2t 2 + 1st) 
+ s 2 m\ ( m| - (s + 2t) mf + 2t(s + t)) 



m. 



(5.3b) 



da j EA,i,j9s{s- ml) \T- 



A\2 



dn 



96 gi g 2 M 2 7r 2 s 5 m^ 



111 \ 3 

( — sj s 



+ 2t ( m~ — s — t) ( ( m~ — s) + 2s ml ml + s 2 m~) 
5.2 Scattering cross-sections 

The cross-section ocm is obtained by integrating the differential cross-section. 
The scattering cross-sections for the E, I and J processes are 



(5.3c) 



2 *rnA\2 

-ji I 



96 0i 02 M 2 vrs 6 m 2 ^ 



+ 2s 2 m~ ( m~ + 3s m~ + s 2 ) 
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(5.4a) 
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36 0i 02 M 2 7rs 5 m 2 ^ 
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S Wlj 
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2( 
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sj m 9 + 4s m^ ( 



m. 



(5.4b) 



72 0! 02 M 2 7ts 5 m^ 



+ s 2 m? ( m? + 2s) 



(m? 



\ 2 4 

s) m g 



+ 2s m? (mj — s) 



m. 



(5.4c) 



6 Gravitino abundances 

Having obtained the cross-section ocm we use the procedure given in section 2 to 
calculate the number density of gravitinos via each process using the integrated Boltz- 
mann equation. We calculate the number density of gravitinos produced via different 
processes E, I and J. For m < T d we only consider t > t t h r - For m > T d , we calculate 
the number density separately from t d to tkin, tkin to t t h r and for t > t thr . Parameters 
needed to evaluate the gravitino abundance are given in section 7 and appendies C 
and D. 
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6.1 m < T d 

As shown in section 4.1 there is no gravitino production until the flat direction con- 
densate decays at tf. The temperature Tf at tf is derived in section 7.7 and equals 
1 x 10 5 GeV . The gravitino abundance generated after tf is given by the standard ex- 
pression for gravitino production in a thermal Universe after inflation with the reheat 
temperature replaced by Tf. For our low value of Tf the gravitino abundance will be 
small enough to satisfy cosmological constraints. 

6.2 m > T d 
6.2.1 td to t kin 

The A terms for the processes E, I, J that enter the RHS of the integrated Boltzmann 
eq. (2.6) for the phase space distribution in eq. (2.4a) are given in appendix B. Below 
we retain only the dominant term for each process on the RHS in the limit that the 
incoming energy m^ (ad/a) is much larger than mg, m g and rrig. 
For the process E 



dnf 
dt 



3Hn d = A E = 72B (m 2 g + 2m?) m\ ( 



72 B (mi 



, 2\ 4 / td 
2m 9> m <t> [ J 



72 x 3Qt _9/2 



(6.1) 



where 



B 



Q 



E 



2 M 2 17^12 p4 
A,i,j iJs m Pl\ 1 ji I 1 d 



256 7r 2 gi g 2 m ~ m , 
Ba^mUlTf 



(6.2) 



(6.3) 



The parameter B is derived in appendix C. We have used also the fact that the gluon 
and gluino mass is time dependent quantities as given by 



m n 



m- g 



a ip 



r„ 



i 



m J \T d t 

as derived in eq. (D.7). Solving eq. (6.1) we then get 

2 



3/4 



216K 



216K 



1- - 



where K is defined as 



K 



Batumi Y d 



2m 2 



t) 

a J 



3/2 



(6.4) 

(6.5a) 
(6.5b) 

(6.6) 
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Similarly for the process I, 



dn~ 



(i) 



dt 



+ 3Hn 5 = A I = 4SBmjml 



\8Bm\ m\ 



48 Qt 



-9/2 



Eq. (6.7) implies that 



48K 



1 - 



(hi 



hi 
t 

ad 



3/2 



(6.7) 

(6.8a) 
(6.8b) 



For the process J, 



96 B m g m q 



96Bm 2 g m\ 



96 Qr 9 / 2 . 



J 
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td 
t 



Eq. (6.9) implies that 



nfit) 



96 K 



9QK 



t) 

Od^ 3 

a 



3/2 



(6.9) 



(6.10a) 



(6.10b) 



One observes that the factor in the square brackets in the above expressions is 
related to gravitino production while the factor ~ 1/ct 3 reflects dilution due to Hubble 
expansion. Because of the fast decline with time of the source term on the RHS of 
the integrated Boltzmann equation, the production of gravitinos shuts off quickly and 
then the number density falls as 1/a 3 . 

6.2.2 tki n to tfar 

For tkin to t t hr we evaluate A defined in eq. (2.6) with the phase space distribution 
function given in eq. (2.4b). 

A = exp(£i/T)exp(£ 2 /T) J dR x dll 2 exp(- J E 1 /T) exp(-E 2 /T)W 12 



exp(£i/r) exp(£ 2 /r) n[ q n eq 
ni(t)n 2 (t){av), 



f tflli dU 2 exp(-E 1 /T)exp(-E 2 /T)W 12 

1 n 2 J 

(6.11) 
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where (av) is the expression in square brackets and 

ni = J exp[-(Si - £i)/T] d 3 Pl for i = 1, 2 



A n! 



Here. 



A % 3 - ( 6 - 12 ) 



Ai = exp(^/T) < 1 for t kin < t < t thr , 

= 1 for t = t thr , (6.13) 



and 



eq 

n- 



J exp{-Ei/T) d 3 Pi , (6.14) 



(2tt) 3 

and we have ignored factors of 0(1). Our expression for rij includes unlike in ref. 
[2]. The number density of all particles in the Universe is 

3 3 

^4. — f7* in thermal equilibrium, where g* is the total number of relativistic degrees 
of freedom for a fully thermalized plasma. For the minimal supersymmetric Standard 
Model (MSSM) = 228.75. Since we can express A for t kin to t thr in terms of the 
standard thermally averaged cross-section we can therefore use expressions for (av) 
from the literature. 

In kinetic equilibrium, temperature is well defined and we choose to re-express 
the integrated Boltzmann eq. (2.6) 

+ 3Hri£ = (av)riin 2 , (6.16) 

as 



IT 

where a sum over all processes E, I, J is implicit. Defining Y' = n/s', where s' 



T-£ + 3Hn 5 = (av) ni n 2 , (6.17) 



2tt£ 4 rp3 
45 



AT , we can rewrite the Boltzmann equation as 



■ dY' 

Since T ~ 1/a, T/T = —H and so 

where we have used H 2 = p/(3M 2 ) and p = (7r 2 /30)*4.T 4 . Expressions for (av) from 
ref. [35], were used in ref. [2]. However these results were for the ±3/2 helicity states of 
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the gravitino. Including the dominant ±1/2 helicity states introduces an enhancement 
factor ~ m|/m~ as discussed earlier. Then 

1 / 777? 

1 / 

<<™>' = 32M? ( X + S| J 16 ° (6 ' 20b) 

2 



/ m 

I 

G 



where we have considered only the SU (3) c gauge group and adjusted the enhancement 
factor for the E process based on eq. (5.2a). Below we shall consider the contribution of 
the ±1/2 helicity states only. The gravitino abundance just before full thermalization 

a t tthr 1S 

y a (i "" ) = 3fer = + £ C Wv)M92MdT (6 - 21) 

where s'{t t hr) and T min are the entropy density and temperature at t t h r just before 
thermal equilibration. s'(t t h r ) = (2tt 2 /A5)AT^ lin and (av) includes a sum over all the 
processes eq. (6.20). Y~(t kin ) is obtained using eqs. (6.5a), (6.8a) and (6.10a) and the 
temperature at t k i n needed for the entropy is given by eq. (7.2). Then 
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in J 



TT'Q (tkin) 
$(tkin) 

Y.lQl K [1 - (td/tkinf] (td/t k 
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(27rV45)^l(m /3) 3 (r fcm /r d ) 3/2 

J2 iqi K[l-(T km /T d ) 2 } 
(2tt 2 /A5)A (m^/3) 3 



(6.22) 



where qi = 216,48,96 for the I = E,I,J processes respectively. Note that K oc B is 
also process dependent. 

After equilibration at tthr the temperature falls to 

T R = - T mm , (6.23) 
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Then, 
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where Y^(t t h r ) = n>a(tthr) / s (Uhr) and s(t t h r ) = (27r 2 /45)g*Tj-j. is the contribution 



/2 s ) 

from gravitino production between tj and tk% n and V- is the gravitino production 

G 

between tkin and tthr- 

In ri 2) , the cross-sections (cry) depend on the enhancement factor m~/m~. The 



G 



gluino mass for this period is obtained in eq. (D.8) and gives 

r\ 2 / pi \ 3/2 / rri 

d \ I J- kin \ ' J 
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, m 0/ V / \Tkin, 

Therefore substituting eqs. (6.20) and (6.25) in eq. (6.24) we have 
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(6.26) 



where /z = 48 x 3/2, 16, 8 from eq. (6.20) for I = E,I, J processes respectively and gi t2 
are as given in appendix C. We have ignored the variation of a with temperature and 
use the value of a = 5 x 10~ 2 . 

6.2.3 t > t thr 

As discussed in section 4.2 gravitino production is suppressed until the flat direction 
condensate decays at £/. The temperature 7/ at tf is given in section 7.7 and equals 
1 x 10 5 GeV . Such a low temperature will lead to a gravitino abundance consistent 
with cosmological constraints. 



7 Obtaining parameters 

The values of parameters which we require to obtain the gravitino abundances are 
evaluated in this section for both mo < Td and mo > Td- 
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7.1 Tkin and T^ in 

Tk in is the temperature at t^n just after kinetic equilibration. Since kinetic equili- 
bration only re-distributes energy among particles without changing their number, the 
average energy per particle is the same before and after kinetic equilibration. Therefore 

\k B T km 

1/2 





a d 
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@>kin 
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(-) 


2 




m 




2 





1/2 

(7-1) 

where the RHS represents the energy of the particles just before kinetic equilibration 
and we presume 2 body decay of the inflaton. The Boltzmann constant ks is hereafter 
set to 1. Therefore 

Tkm = ? (t?) ' (7 - 2) 

where Tkin is given in eqs. (4.5b) and (4.21). 
For m < T d , 

-\aM 2 T 1/2 m 3/2 
^ n =* a Pl \ ° =UGeV (7.3) 



For m > T d , 



kin 



m 4> 



3 a Mj, t ml 



2 m^ipl 

= 1 GeV (7.4) 

Thus, the value of T^ n for mo < is 1 x 10 11 GeV and for mo > is 1 x 10 12 GeV . 
7 9 T 

' - 1 mm 

The temperature T min just before thermalization at t tnr can be obtained by 

akin 



± rain - 1 - kin 



J kin 



O-thr 

tkin 
tthr 



1 

( ^ ) 2 (7-5) 



2 x 10 10 GeV for m < Y d 
2 x 10 11 GeV for m > Y d 

where we have used the relation T t hr = o^kin as discussed in section 4. 
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7.3 A 

To obtain A we note that the energy density of the Universe Pki n + just after kinetic 
equilibration at tki n is equal to the energy density pkin- just before kinetic equilibration. 
Then 



or, 




UxlO- 5 , for m <T d = 10 4 GeV 
or, ^4 ~ < „ (7.7) 

[4 x 10" 11 , for m > T d = 10 GeV. 

Thus A<^ 1. 

7.4 A 

The number density of a species i just after kinetic equilibration is Aig{Tl in l -a 2 . This 
equals the number density of the species just before kinetic equilibration which, in 
turn, is the number density of the species at t d scaled by 1/a 3 . The decay rate of the 
inflaton to a species i will be proportional to gi and so ni(t d ) is also proportional to 
gi. Then At is the same for all species. Since A = J2j ^j9j = Ai Ylj 9j ~ 200 Ai and 
At « .A/200. 

7.5 T« 

As stated earlier in eq. (6.23) the temperature at t thr is given by 

t - ^V /4 t 

-L R I I -L min 

\9*J 
[A\ 1/A 

= ( — j V 7 " Tjfcin 
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For m < T d , 



T, 



R 



\/6 W ^ ry 4 

5xl0 8 GeV. (7.9) 



1 f A\ l/i am]l 2 M Pl 

Tr 



For m > r d , 



= 2xl0 8 GeV. (7.10) 

where we have used eqs. (7.3) and (7.4) for the expression of T kin for m < T d and 
m > T d respectively, and a = 5 x 1CT 2 . 

7.6 T G 

The temperature Tq at i ff is given by 

Ta = T R (^Y (7.11) 

For m < r d , using eq. (4.13) to define og, we get 

T 3 

7b — 2 ( m otthr) 



^\ 3 / 4 m 3 m 3/2 



""o 



27 ;W a^rf 
= 300GeV. (7.12) 

where we have used section 4.1.2 for t t h. r and eq. (7.8) or (7.9) for Tr and a = 0.1. 
For m > T d , using eq. (4.24), the temperature at to is 



1 / ,4\ 3/4 



27W «^o^ 
= 1000GeV. (7.13) 

where we have used section 4.2.1 for t t h. r and eq. (7.8) or (7.10) for Tr and a = 0.1. 
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7.7 T f 

The temperature at tf is given by 



Tf - T R (t thr /t f ) 



1/2 



(7.14) 



For m < T d , 



1 = i mv 7 v 

7 3 \gj J 




= 1 x 10 5 GeV . 



(7.15) 



where the expressions for t t h r , tf and Tr are obtained from section 4.1.2, eqs. (4.16) 
and (7.8) or (7.9) respectively. 
For m > T d , 



where the expressions for tthr, tf and Tr are obtained from section 4.2.1, eqs. (4.26) 
and (7.8) or (7.10) respectively. 

8 Results 

For the case m < T d , with m = 100 GeV, T d = 10 4 GeV and tp = 10 14 GeV, 
t d , t ,t kin , t thT and tf are lO^GeV" 1 , lO^GeV" 1 , 9 x lO^GeV" 1 , 2 GeV- 1 and 

4 x 10 7 GeV _1 respectively, tp is constant at 10 14 GeV from t d to £ an d falls to 
2 x 10 13 GeV and 2 x 10 12 GeV at tki n and tthr respectively. The gluon and gluino 
mass is constant from t d to t at 2 x 10 13 GeV and falls to 4 x 10 12 and 4 x 10 11 GeV 
at tkin and t t hr respectively. The final reheat temperature Tr after full equilibration is 

5 x 10 8 GeV, which is less than the gluon and gluino mass at that time. As we have 
argued in section 4.1 there is no gravitino production in the non-thermal phase, nor 
after thermalization till the flat direction condensate decays at tf, because of phase 
space suppression associated with the large mass of the outgoing gluon or gluino. The 
temperature Tf is 1 x 10 5 GeV and gravitino production after tf is hence small and 
within cosmological constraints. 

For m > T d , with m = 100 GeV, T d = 10 GeV and (p = 10 14 GeV, t , t d , t fein 
and t thr are 0.01 GeV" 1 , 0.1 GeV" 1 , 0.9 GeV" 1 and 20 GeV' 1 . <p faUs to 10 13 GeV at 
t d , 2 x 10 12 GeV at tkin and 2 x 10 11 GeV at t t hr- The corresponding gluon and gluino 
masses are 2 x 10 12 GeV , 4 x 10 11 GeV and 4 x 10 10 GeV at t d , t Hn and t thr respectively. 
The final reheat temperature Tr is 2 x 10 s GeV . The gravitino abundances generated 
in the different epochs are 




1 x 10 5 GeV . 



(7.16) 



Vi 1} = 3 x 10 3 

G 

Yi 2) = 1 x 10' 



-2 



(8.1) 
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where Y~ is the abundance due to gravitino production between td and t^ in and 

(2) 

Y~ is the abundance due to production between tkin and t t h. r - Gravitino production 
after thermalization at tthr is suppressed until the condensate decays at tf. Since the 
temperature at tf is 1 x 10 5 GeV the gravitino abundance generated after tf will be 
small and within cosmological constraints. 

The gravitino abundance for mo > (in each epoch before t t h r ) is larger than the 
upper bound on the gravitino abundance from various cosmological constraints. The 
source of the large abundance is the enhancement by the factor ~ m~ (</?)/ m~ which 
was overlooked in earlier works. Therefore, this scenario with m > 1^ is ruled out 
due to the over- abundance of gravitinos for the parameters above. 

We have chosen tpo = 10 14 GeV such that gravitino production commences from 
td itself so as to test whether a non-thermal Universe can lead to a low gravitino 
abundance. However we have obtained a large gravitino abundance for mo > Yd due 
to the large gluino mass. But one can also suppress gravitino production using a 
large gluon and gluino mass. If we set t p in eq. (4.27) greater than t thr , gravitino 
production is phase space suppressed till t t h r - We find that this condition is satisfied 
for (po > 3 x 10 15 GeV . Furthermore, eqs. (4.25) and (4.26) will imply that gravitino 
production after t t hr will also be disallowed till the condensate decays for such large 
values of (p . 

9 Conclusions 

In this article we have studied the novel scenario proposed by refs. [1, 2] where the vac- 
uum expectation values of SUSY flat directions generate a SUSY conserving large mass 
for gauge bosons and gauginos. Due to the large gauge boson masses thermalization 
is delayed after the end of inflation. In the non-thermal Universe gravitino production 
in refs. [1, 2] is suppressed because of the dilute nature of the plasma. The delay 
in thermalization also leads to a low reheat temperature and so gravitino production 
after thermalization is also suppressed. Thus the gravtino problems of SUSY models 
is avoided in this scenario. 

In this article we have used the relevant Boltzmann equation for the non-thermal 
scenario and have obtained the gravitino abundance due to processes involving incom- 
ing quarks and squarks and an outgoing gravitino and a gluon or a gluino. We have 
investigated two cases (i) where the flat direction starts oscillating after inflaton de- 
cay (m < Td) and (ii) where the flat direction starts oscillating before inflaton decay 
(m > T d ). 

In the first case (m < Td) we find that there is no gravitino production during 
the non-thermal phase because of phase space suppression due to the large gluon and 
gluino mass. Gravitino production is suppressed even after thermalization till the flat 
direction condensate decays at tf. However the temperature Tf is low and as such any 
gravitino abundance produced after tf is small and is within the standard cosmological 
bounds. 

In the second case (mo > Td) there is large gravitino production before thermal- 
ization. This is because the ratio m~/m~ which appears in the (helicity ±|) gravitino 
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production rate is large. Therefore gravitino production is greatly enhanced and leads 
to a conflict with observational constraints on the gravitino abundance. (However 
the gravitino abundance can be decreased by increasing the initial value of the flat 
direction vev, ip , leading to phase space suppression of gravitino production.) After 
thermalization, as in the case m < T^, there is no gravitino production before the flat 
direction condensate decays. After the decay of the condensate the temperature Tf is 
low and does not lead to further conflict with cosmological bounds. 
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A Derivation of the phase space distribution functions 

The number of incoming particles Ni for gravitino production is given by 



For the time interval t d < t < t kin the form of / is given by eq. (2.4a). Substituting 
that in the above equation we get, 



where we have considered E{ ^> rrii. Now at t = td, the number density of the 
particles is equal to the number density of the inflaton, as — > Xi + X 2 . So, riiit d) = 
n<t>(td) = P4>(td)/ m (j>- Then for td < t < t kin , the number density riiit) is given by 
riiit) = (p^itd) lm^){a d j a) 3 . Substituting the above expression for rij(t) in eq. (A. 2) 
we get 



B A terms for td < t < tkin 

Below we present the A terms for the processes E, I, J during the interval t d < t < t kin 
as needed for the m > Y case using eqs. (2.6), (2.8) and (2.9), and the phase space 
distribution in eq. (2.4a), with the expressions for <Jcm as given in eqs. (5.4a), (5.4b) 
and (5.4c). 




(A.l) 




(A.2) 




(A.3) 



29 



The A term for E process is 



A E (t) 



6x2x2xE AM CiC2g s 2 |^| 2 
18432a 4 M 2 7i 5 a 8 d mlml 



[3 m g m^a 12 + 4 a 2 , m 6 g m| (3 rr 



9 9 
ml 



Ami) 



10 m~ m g 



+ 18 a d m A g m\ ( m A g - 4 m\ m 2 g + 2 m|) mja 8 + 12 a° m 2 , ( 
+ 18m|m 2 -4m?) m^a 6 + (-9 mj + 36 m? m* - 120 m\ m 2 g + 25 m 



6 In 



a 2 m 2 



(3 mjj — 24 m 2 ml + 24 rr 



36 ai°mj( 



2 m 2 ,) r 



9 9 

10„2 



2 ml) 



nr.a 



3a 12 (m 2 g + 2ml) mf] 



(B.l) 



where T A is the generator of the SU(N) group in the fundamental representation 
and Ci oc (ad/ a) are given in eq. (A. 3). A E includes a factor of 6 for 6 flavors, a 
factor of 2 for the 2 hermitian conjugate E processes and another factor of 2 for the 2 
squark 'chiralities'. Retaining the dominant term in the above expression in the limit 
m^dd/a) 3> mg, m g and m g we get 



A E (t) = 72B (m 2 + 2m 2 ) 



rrij 



where the parameter B is as defined in appendix C. 
The A term for the I process is 



Mt) 



6912 M 2 n 5 m% 



(40 



m. 



36 m~ m 2 



9m|) m~ 



2a 6 



(B.2) 



4 8 

m g m~ 



+ 



+ 



2a 4 m 2 (2 



m; 



2a 2 (6 m 4 — 6 m~ m 2 + m~) m 



4 ™4 ^,2 



4 4 



a 2 , m 2 



2a 4 



m- 



2 a\ m 2 ( m 2 — 2 m~) m; 



(B.3) 

A/ includes a factor of 6 for the flavors. The term which is highest order in m^a^/a) 
in the above expression is 

,fi ' 



Mt) = 

The A term for the J process is 



485 m\ m\ 



(B.4) 



6x2xE,, y C 1 (7 25s 2 |^| 2 



6912 a 4 M 2 vr 5 a 6 , ml 



m G m 4 



m 9 m^a + 3 a d m fl m^ y m g — 2 m g ) m^a 



+ 3 a d m| (6 m A - 6 m~ m 2 + m|) 



mia 6 



+ a d mj 



- 10 mt 



9 m^ m g — 9 m^ + 6 In 



a^m: 
a 2 m 2 g 



3^( 



m„ 



2 m? m 2 



3-1) 



m 8 ,a 2 



3 m| m 



in 



(2 m 4 - 3 m 2 m?) ] m\a A 
(B.5) 



- 30 - 



A j gets a factor of 6 for flavors and a factor of 2 for the squark 'chiralities'. Similarly, 
as before, the term highest order in m^(ad/a) for the J process is 



Aj(t) = 965m?m 



(B.6) 



C B parameter 

To simplify our equations we define new parameters called B and B for the pre-factor 
present in A terms, i.e., 



B 
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18432 M 2 pl 7T 5 m| 



(c.i; 



where M Pl = 87rM 2 . Substituting the values of C\ and C 2 from eq. (A. 3) in the above 
equation gives us 



B 



P*{td) (2tt) 



7l 



9i92 



E 



j 1 1 



18432 M 2 pl vr 5 m| 



(hi 



B 



"a 



(C.2) 



where 



( 27r ) 



7T 



0102 



18432 M Pl vr 5 ml 



G . 



m' 



7T 



8tt 



A|2 



V o 2 IT 4 

/—/A,i,j its \- L ji 

18432 5l g 2 7T 5 m| 



256 7T 2 5-i 02 ™~ ™J 



(C.3) 



where we have used p^{t d ) = (3/ (8n))M Pl Y 2 d . We suppress the index E,I,J for 
associated with the different g 12 for each process. Considering the spin and color 
of the incoming particles, (01,02) = (6, 3), (6, 6) and (3,3) for the E,I,J processes 
respectively. 



D Mass of gluon and gluino 

The mass of the gluon and gluino m g j ~ which varies as 1/ a 3 / 2 after the condensate 
starts oscillating at t . Below we derive the mass of the gluon and gluino for various 
epochs in time for the two cases, m < and m > IV 
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D.l m < T d or t > t d 

1. For td < t < t , ip = (p . So 



m 9j g{t d ) = y/aifQ . (D.l) 



2. For t <t < tthr, 



m 9, 9 = "igM (jf) 



y/aip ( j 



3/2 
3/4 



\ 3/4 



&*[— t ) (D.2) 



At t — tki n , 



m g ~ g {t kin ) = y/a(p ( , (D.3) 

1 mo 



where T^n is given by eq. (7.3). 

At t = tth r , 

1 tftr 



m g ~ g {t thr ) = y/aipo ( — - ) • (D-4) 
m 



3. For t > t t hr, 

'ath. r \ 3 / 2 



m g , g = "WW) ( — J 

r fr thr \ 3/i (T\ i/2 



D.2 7ii > T d or t < U 
1. At t = t d , 



m 



</■</ 



g(td) = 



V^O ( — ) 

\adj 



3/2 



= yfatpo (E±) , (D.6) 

where we have used the fact that the Universe was matter dominated for the 
period to < t < td and a ~ t 2 / 3 , as relevant for an inflaton oscillating in a 
quadratic potential after inflation. 
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2. For t d < t < t ki 




(D.7) 



where we have used the fact that the Universe is radiation dominated for the 
period td <t < t^n an d a ~ t 1 / 2 . 

3. For thin < t < tthr 




where we have used the relation a ~ 1/T for t^n <t< tthr- 
4. For t > t thr 




(D.9) 



where we have used eqs. (D.8) and (7.5), and a ~ 1/T for t > t thr . 
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